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Dynamical dark energy (DE) phenomenon emerges as a geometrical effect accompanying the 
cosmological expansion of nonrelativistic fermionic matter. This occurs without the need for any 
fluid, like e.g. dynamical scalar field (quintessence, cosmon, etc.), and with conventional form of 
the Einstein equations in contrast to other known geometrical DE models. The phenomenon results 
from first principles in the framework of the two measures field theory where, in the Einstein frame, 
both fermion masses and the cosmological constant turn into functions of the cold fermion density 
n. This n dependence becomes negligible in regular (laboratory) conditions but it has a crucial role 
in cosmology: the cosmological dynamics is governed by the scale factor dependent fermion density 
^S) , yielding the results similar to what quintessence does. In the 4D gravity model where the original 

T-H ■ action involves only the cosmological constant and massive fermions without selfinteraction, in the 

■ broad range of the parameter space, we find a rapid transition (at z ~ 1) from the cold fermionic 
04 , matter dominated epoch with the equation of state (EoS) w = p/ p « to the accelerated epoch 

^jrj with w < —1/3 that is followed by a further asymptotic decrease towards w — —1. We argue that 

■ the fermions should be associated with neutrinos. Depending on the parameters, the late time 
I evolution can proceed either with monotonic decreasing EoS to w = —1 or in the phantom- like 
. regime. In the latter case, for a broad range of initial conditions, w crosses the phantom divide due 

\^ ■ to a new type of the neutrino DE effect where neutrinos pass through the state with zero mass and 

^S) ' with the vacuum-like EoS — —p,j. After passing their minima, w and p increase asymptotically 

, with cosmic time approaching from below w = —1 and p — const respectively and thus providing a 

pseudo-rip scenario. 

[ PACS numbers: 04.20.Cv, 04.50.Kd, 95.36.+X 

^ I I. INTRODUCTION 

>: 

04 . Numerous attempts to understand the nature of the discovered cosmic acceleration [Tj have originated the common 

■ feehng in the cosmology and particle physics community that even though all observational data are in agreement with 

■ the simplest model where this phenomena is described by the cosmological constant (CC), the latter is not enough to 
give adequate answers to all challenges of the cosmology. Therefore to explain the origin of the acceleration it seems 

, to be impossible to advance without some kind of changes in fundamental of modern physics 0. The new physical 

■ entity responsible for the present accelerated expansion of the universe was denominated as the dark energy (DE) 
I [ and the problem was exposed intensive study with vast of models, for reviews seeQ-Q. Using terminology of the 

• • ■ reviewQ, to a large extent, all models with variable DE can be divided into two main classes: 1) physical DE models 

. ^ ! where DE is associated with some fluid usually realized by means of a new, very weakly interacting physical field. 2) 

K> . geometrical DE models where effect of DE is realized by means of modified gravity. 

^ ■ In the most of the physical DE models, like e.g. quintessence 0, coupled quintessence mass varying neutrinos 



in tne most or tne ptiysicai ut^ moaeis, uKe e.g. qumtessencepj], coupiea qumtessence[8||, mass varymg neutrmos 
and neutrino DE Q-fll], growing neutrinos [isj. chameleon cosmology [l^. phantom [T5|. quintom[l6[, the DE entity 
is modelled by scalar field(s) with very special features (very flat potential, etc.) or/and using non-scalar systems, 
e.g. models with specially adapted spinors[l7[ . Using of such fields means the exit from the framework of the matter 
content of the particle Standard Model. 

The geometrical DE models, e.g. higher dimensional DGP braneworld DE [lB|,[3:^CR) DE scalar-tensor 

,11, P 



mo dels [2 1|. braneworld with Einstein-Gauss-Bonnet gravity [22j. etc. (see also reviews [3|,[5|,[6|) demonstrate possi- 
bilities to realize an effective DE without any field or other matter entity intended to produce DE. However the price 
for this effect is the emergence of nonconventional 4D gravitational dynamics. Only after certain regrouping of terms 
having a geometrical meaning to the effective matter source terms, the gravity equations takes the Einsteinian GR 
form. 

The model of the present paper is realized in the framework of the two measures field theory (TMT) which is a 



kind of alternative theories [23[. In this respect we want to stress that TMT possesses a number of peculiar features 
which allow to regard it as a new type of alternative theories. First, the modification of the gravity and matter 
sectors in TMT is implemented on the common basis. Second, apart from proceeding in the first order formalism, 
all the modification in the action of TMT , as a matter of principle, reduces to the use of the volume element with 
two measures: the regular one ^/—g and the new one $ involving additional non-dynamical degrees of freedom (for 
details see Sec. II). Third, TMT admits but does not require to involve in the underlying Lagrangian terms different 
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from those of the conventional theory (i.e. Einstein gravity and Standard particle model). Fourth, in the Einstein 
(physical) frame, the gravity and all matter field equations have conventional form without the need in regrouping and 
renaming geometrical and matter source terms. The novelty consists in emergence of variable cosmological constant, 
very nontrivial form of scalar field potential, fermion selfinteraction and Yukawa type coupling constants become 
local function of the fermion density. Therefore fermions generically may possess very unconventional properties that 
makes it possible, for example, to realize effect of the neutrino driven DE which is the subject of this paper. But 
properties of fermionic matter in regular conditions are undistinguished from those in the conventional theory. The 
term 'regular conditions' means here that the local fermion energy density p/ is tens orders of magnitude larger than 
the vacuum energy density Pvac, that is actually a condition which is fulfilled in all measurements carried out so far in 
physics, including tests of GR. For example in regular conditions, corrections to fermion mass have the order p^ac/ Pf- 
In this respect, in contrast to other alternative theories, the novelty of the most of the studied so far TMT models is 
that in order to achieve for undesirable deviations from the conventional theory to be unobservable there is no need 
in special tuning or constraints on the masses and coupling constants in the underlying Lagrangian. 

In our earlier publications on TMT[23l-[23j[i3 we concentrated on studying the scale invariant TMT models. Here 
the scale invariance is spontaneously broken by the integration of the equations of motion and we find that the dilaton 
plays the role of the inflaton field in the early universe and that of the DE field in the late universe. It was shown [28| 
that the Yukawa type effective coupling constant / of the dilaton to matter depends on the local matter density. For 
matter in regular conditions the magnitude of / is suppressed by the factor of the order of Pvac/ Pf- The model yields 
such resolution of the fifth force problem automatically, without any special tuning of the parameters. 

In Ref. [13] devoted to cosmological applications of the model, we have studied the effect of neutrino DE emerging in 
the course of evolution of the late time universe filled with the homogeneous dilaton field and the cold gas of uniformly 
distributed non-relativistic neutrinos. A new kind of regime in the cosmological dynamics have been revealed where 
neutrinos undergo transition to a state with unbounded grows of their mass and this process is accompanied with 
reconstruction of the scalar field potential in such a way that this dynamical regime appears to be energetically more 
preferable than it would be in the case of the universe with no fermions at all. 

Recently another TMT model [s^ involving a gravitating scalar field with a non-conventional kinetic term as in the 
string theory tachyon and with an arbitrary potential have been constructed. This model gives a unified picture of 
dark energy and dark matter. 

In the present paper we work with the TMT model which may be regarded as a simplified TMT modification 
of the Standard Model. This is the case because even though the underlying Lagrangian of our simplified model 
involves only free massive fermions and the CC term, the TMT allows [26j to extend it by adding all fields and 
symmetries of Standard Model without changing the main results of the paper. The reason is that the key role in 
generation of the aforementioned effects of the novelty of TMT (in comparison with conventional theory) belongs 
to the algebraic constraint which appears due to the new measure of integration $. The constraint determines the 
scalar field C, — ^/^/—g as a local function of \E'\I/ (where ^I' refers to massive fermion fields) and the latter is the only 
contribution to the constraint from the Standard Model fields. Therefore, at least in the first step of studying the 
cosmological applications of the model, one can restrict ourself with the simplified model. 

In spite of the fact that aside from fermions, our simplified model does not contain any additional dynamical degree 
of freedom, the variable effective CC Atot is generated and its local space-time value is determined by the local fermion 
density. In this paper we explore the cosmological dynamics of the late spatially flat FLRW universe starting from the 
cold matter domination epoch. Atot is governed by the density of the cold fermion gas involved in the cosmological 
expansion. Transition to the Atot domination epoch happens very rapid (about z ~ 1) and this is achieved without 
any special tuning of the Lagrangian parameters. Apart from the evident possibility of the late time scenario with a 
monotonic decreasing of the equation of state (EoS) w — > — 1, the model admits also the phantom-like regime. The 
main goal of this paper is to demonstrate that the Standard Model formulated in the framework of TMT generates 
variable CC without the need for additional dynamical degrees of freedom. To our knowledge, this is the flrst model 
where this idea is realized, see for example Refs.fsij.fs^. 

The organization of the paper is the following. In Sec. II we present: a review of the basic ideas of TMT; discuss 
properties of the new degrees of freedom; describe on qualitative level some details of using the action principle 
resulting in the equations of motion and the constraint in terms of the original set of variables; explain, again 
qualitatively, the reasons and the results of transition to the Einstein frame. This allows us, after the formulation 
of the model in Sec. Ill, to present the equations of motion in the Einstein frame at once. Sec. IV is devoted to the 
very detailed analysis of physically permitted intervals for the geometrical scalar field C- In Sec.V we describe new 
aspects of the cosmological averaging procedure caused by the appearance of the constraint. Afterwards we present 
equations of cosmological dynamics resulting from our field model after cosmological averaging and classify possible 
scenarios of the late universe. In Sec. VI the results of numerical solutions are collected and classified in accordance 
with qualitative analysis of Sec.V. In Sec. VII we shortly summarize main results of the model that should help to the 
reader to see the picture altogether and finally present some preliminary speculations concerning the possibility to 
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describe the fermion dark matter as a distinctive fermion state in TMT. 



II. MAIN IDEAS OF THE TWO MEASURES FIELD THEORY 



A. Two volume measures 

TMT is a generally coordinate invariant theory where the action may be written in the form(23j 

J Li$dlx + J L2V^d'^x (1) 

including two Lagrangians Li and L2 and two measures of integration: the usual one ^/—g and the new one <f> 
independent of the metric. The measure $ being a scalar density may be defined as following: 

^ = e^"'°''^eabcddf,ipadu'fbda'PcdfHPd, (2) 

where ipa{x) (a = 1, .., 4) are scalar fields. Since $ is a total derivative, the shift of Li by a constant, Li — > Li +const. 
has no effect on the equations of motion. Similar shift of L2 would lead to the change of the constant part of the 
Lagrangian coupled to the volume element ^J—gd^x. In standard GR, this constant term is the cosmological constant. 
However in TMT the relation between the constant term of L2 and the physical cosmological constant is very non 
trivial. 

In addition to the above idea concerning the general structure of the action in TMT, there are only two basic 
assumptions: 

(1) The Lagrangian densities Li and L2 may be functions of the matter fields, the metric, the connection (or spin- 
connection ) but not of the "measure fields" {pa- In such a case, i.e. when the measure fields Lpa enter in the theory 
only via the measure the action ([Ij has the infinite dimensional svmmetrv(23|: ^Pa ^ + fa{Li), where fa{Li) 
are arbitrary functions of Li. 

(2) One should proceed in the first order formalism where all fields, including metric and connections ( or vierbeins 
ej^ and spin-connection ujf in the presence of fermions) as well as the measure fields (pa are independent dynamical 
variables. All the relations between them follow from equations of motion. The independence of the metric and the 
connection in the action means that we proceed in the first order formalism and the relation between connection and 
metric is not necessarily according to Riemannian geometry. 

Varying the measure fields (^q, we get 

Biid.L, = (3) 

where 

- e^-^^^eabcdd^Vbdc^^cd^Vd. (4) 
Since DetiBi^) = it follows that if $ 7^ 0, 

Li = sM'^ = const (5) 
where s = ±1 and M is a constant of integration with the dimension of mass. 



B. Generic features of TMT 



First of all one should notice the very important differences of TMT from scalar-tensor theories with nonminimal 
coupling: 

• In general, the Lagrangian density Li (coupled to the measure $) may contain not only the scalar curvature 
term (or more general gravity term) but also all possible matter fields terms. This means that TMT modifies 
in general both the gravitational sector and the matter sector in the same fashion] 

• If the field $ were the fundamental (non composite) one then the variation of $ would result in the equation 
Li ~ Q instead of ([5]), and therefore the dimensionful parameter M would not appear. 
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Applying the first order formalism one can show (see for example Ref.[2^) that the resulting relation between metric 
and connection, as well as fermion and scalar fields equations involve the gradient of the ratio of the two measures 

(6) 



which is a scalar field. For understanding the structure of TMT it is important to note that TMT (where, as we 
supposed, the scalar fields enter only via the measure $) is a constrained dynamical system. In fact, the volume 
measure $ depends only upon the first derivatives of and this dependence is linear. The fields do not have 
their own dynamical equations: they are auxiliary fields. All of their dynamical effect, at least at the classical level, 
is displayed only in the following two ways: (a) in generating the constraint; (b) in the appearance of the scalar field 
C, and its gradient in all equations of motion. 

The origin of the constraint is clear enough. There are two equations containing the scalar curvature: the first 
one is Eq.(l5]) and the second one follows from gravitational equations. The constraint is nothing but the consistency 
condition of these two equations: eliminating the scalar curvature from these two equations one obtains Q{x) as a 
function of matter fields. The surprising feature of the theory is that neither Newton constant nor curvature appear 
in this constraint which means that the geometrical scalar field C,{x) is determined by the matter fields configuration 
locally and straightforward (that is without gravitational interaction). 

The appropriate change of the dynamical variables, which includes a transformation of the metric, results in the 
formulation of the theory in a Riemannian (or Riemann-Cartan) space-time. The corresponding conformal frame 
we call "the Einstein frame". The big advantage of TMT is that in the very wide class of models, the gravity and 
all matter fields equations of motion take conventional GR form in the Einstein frame. All the novelty of TMT in 
the Einstein frame as compared with the standard GR is revealed only in an unusual structure of the cosmological 
constant, scalar fields effective potential, masses of fermions and their interactions with scalar fields as well as in the 
unusual structure of fermion contributions to the energy-momentum tensor: they appear to be C, dependent. This is 
why the scalar field C,{x) determined by the constraint, has a key role in these effects. 

III. THE FIELD THEORY MODEL 

Our simpificd TMT model deals with the selfconsistent system involving 4D gravity, the cosmological constant 
(CC) and fermions. Keeping the general structure of Eq.([T]) it is convenient to represent the action in the following 
form: 



b-S 



i j d'^x{^ + by^)R{uj, e) - y d*xy/^AQ 

i 

- J d4x($ + hV^) (7) 

■i 

where ^'^ are the fermion fields of species labeled by the index i; fXi are the mass parameters; = d + ^uj^dcd (where 
we ignore gauge fields); k — 167rG; R{uj, V) — e°-^^e'^'' R^^^abi^) is the scalar curvature; and w^'' are the vierbein and 

spin-connection; g^^'^ = e^^e'j^r]"-'' and Ri_iuab{ij-j) — dmut^i^ab + oj^a'^vcb — (^i O i'); we use the signature (-1 ); Aq is 

a constant of the dimensionality (mass)'*' which in the standard GR would be the CC. Notice that there is no need 
for including similar term with the measure $ because such a term is the total derivative and does not contribute 
into equations of motion. Constants b,k,h are non specified dimensionless real parameters and we will only assume 
that b > and they are different but have the same order of magnitude 

br^kr^h. (8) 

We would like to stress that except for the modification of the general structure of the action according to the 
basic assumptions of TMT, we do not introduce into the action ([7]) any exotic terms and fields. Without changing 
the results of the present paper, one can generalize the model to non-Abelian symmetry adding also gauge and Higgs 
fields to reproduce the standard model (see [1^). However, for purposes of the present paper it is enough to restrict 
ourself with the simplified model involving only massive fermions and the CC term in the action. We would like to 
lay emphasis that refusal from a dynamical scalar field, like quintessence, means that our model does not exit the 
matter content of the standard model. 
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Variation of the measure fields Lpa yields Eq.([5]) where ii is now defined, according to Eq.(IT]), as the part of the 
integrand of the action JT]) coupled to the measure $. 

Equations of motion for fermion fields as well as the equations for the spin-connection resulting from ([7]) in the 
first order formalism contain terms proportional to d^C, that makes the space-time non-Riemannian and equations of 
motion - non canonical. However, with the new set of variables 



(C + hf/i 



(9) 



which we call the Einstein frame, the spin-connections become those of the Einstein-Cartan space-time with the metric 
5^1/ The fermion equations also take the standard form of the Dirac equation in the Einstein-Cartan space-time where 
now the fermion masses become C, dependent 



m,(C) = ii.,F{C), where F{C) 



iC + h) 



We have supposed here that 



C + 6> 



and 



C + fc > 0, 



which will be assumed in what follows. 

It is easy to check that the gravitational equations in the Einstein frame take the standard GR form 



(10) 



(11) 



(12) 



where Gfj,u{gai3) is the Einstein tensor in the Riemannian space-time with the metric gfj,i^] the energy- momentum 
tensor T^;^^ is now 



7^e/7 _ rp{J.can) , rp{A) 



(13) 



where j'^^''^"") ig the canonical energy momentum tensor for fermions in curved space-timefssl]; T^^}^ is the total 
variable CC term in the presence of fermions 



where 



and 



Tj^t^ = ~9^,uAtot ; Atot = Kf f (C) - A,^^^ 



dyn'' 



Ae//(C) 



bsM^ - Aq 



Al;U^(c)^^(c)E/^^*"^^ 



(14) 



(15) 



(16) 



z{0 



(C-Ci)(C-C2) 
2(c + fc)(C + M 



(17) 



Cl.2 - I ^ 



k-3hi± v/(fc - 3/i0^ + 8&(fc - h,) - Akh, . (18) 
Taking into account our basic assumption ([5]) concerning the orders of the parameters 6, k and h we see that 

Ci^C2--b^k^h. (19) 
if no a special fine tuning of the parameters is assumed. 
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The first new effect we observe here consists in the appearance of the variable CC Aeff{C) generated by the 
integration constant sM* (see Eq.®) and by the Aq term of the original action ([7]). Ae//(C) is positive provided 

bsM^ > Aq. (20) 

The second novelty consists in generating the fermion noncanonical contribution into the energy momentum tensor 

rp{f,noncan) _ _~ A (/) /r,i \ 

which has the transformation properties of a CC term but it is proportional to ^'^VE'^. This is why we refered to itflT 
as "dynamica 
pressure and 



as "dynamical fermionic A term". The appearance of means that even cold fermions generically possess 



pinoncan) _ _ (noncan) _ a (/) _ (C Cl)(C C2) V"^ ^ , 

where jg i^j^g noncanonical contribution of fermions into the energy density. 

The scalar field C is determined as a function of ^'^^'^ by the following constraint 

(6 - C)sM^ - 2Ao _ inoncan) 

In the presence of fermions, the r.h.s. of the constraint is nonzero, and therefore C. 7^ const. One can then expect 
that the variable CC term Tjit^ ,Eq.^, mimics the dynamical DE effect in the late time universe. This is indeed 
the case as we will see later on in the paper. Using the constraint (j23p . T^^^* may be represented in the form 

T^t^ = g,.AUO, AtotiO = (24) 

where only the ( dependence remains explicitly. 

Let us imagine for the moment that we were succeeded to solve analytically the constraint (I23p . which is the 5-th 
degree algebraic equation with respect to (C + 6)^/^. Having the function C = C {J2i Mi^'i^O ^^id inserting it to the 
fermion mass mi{(^), Eq. ljlOp . we would obtain a local fermion (very nonlinear) selfinteraction in the Dirac equation. 

Inserting the same solution for C. into Tl^t\ Eq. ([M)) . we would obtain that the DE effect in the late time universe is 
governed only by the expession X^j/^i^'i^iJ while the latter is locally determined by the dynamics of fermions with 
the described selfinteraction. In this respect, the key point of the TMT model under consideration consists in the 
absence of a DE fluid. And even though Atot has the form of the DE-like effective potential it is not a potential of 
any field. 

The described approach runs into the mathematical problem when we want to solve analytically the 5-th degree 
algebraic equation. Instead, one can proceed with the scalar field ^ evaluated by means of numerical solving the 
constraint. To realize the program of numerical solution, it is useful to notice that in what follows we are going to 
study the application of this field theory model to the case of nonrelativistic fermions in the cold fermions (dust) 
approximation. This means that we neglect the effect of fermion 3-momenta. The only component of the canonical 
fermion energy-momentum tensor y^^'"^""^ which survive in this approximation is the energy density Xqq'^'^"'^ — p"j°'^^\ 
Making use of the Dirac equation in the same approximation we obtain 

p(-)=^™,(C)*^.*^ (25) 

i 

To provide the consistent description of the fermion energy density and pressure we have to regard ^'i^'i as the 
field operator. Then the standard semiclassical approach to the Einstein Eqs.([T2]) (when the gravity is treated as the 
classical field) implies the need to evaluate the matrix element of ^''i^'^ between appropriate states of the Fock space 
resulting in the fermion (particles and antiparticles) number densities and then 

)=^„i,(C)n. (26) 
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In this respect, the only new element in our model is the appearance of the geometrical scalar field ^ (treated in 
semiclassical approach also as the classical field) and the constraint . Proceeding with the constraint in the same 
way as with the Einstein equation, we obtain 

{b-OsM^_2Ao ^ pinoncan) ^ where P)"") = A^-J,!^ = Z(C)p^'™^ (27) 

In the case of nonrelativistic fermions, corrections to the cold fermion approximation reduce to the appearance of 
the Fermi momentum depending factors 11 1, jsj] in front of rii in Eqs. (PB|) . (P7| . For the goals of the present paper 



it is enough to notice that these factors are positive. 

Variability of the CC term ([M)) is governed by the geometrical, nondynamical scalar field Ci^)- In its turn, due to the 
constraint, the local space-time value Cix) is determined straightforwardly (without gravity) by the local magnitude 
of the fermion density . And vice versa, the dynamics of the fermion field is affected (again locally) by C(x) via the C 
dependent fermion mass, Eq. (fTO| . 

This way of solution to the problem does not affect the above conclusion about the absence of a DE fluid but it 
is rather technical method to get round the pure mathematical problem. In other words, the emergence of ^ in the 
numerical solution to the problem does not mean that the model contains a scalar field describing DE. 

IV. PHYSICALLY PERMITTED INTERVALS OF ( 

In the TMT models every local magnitude of the fermionic matter density determines the local value of the ratio of 
two measures C = ^/^—g through the algebraic constraint (j23p which does not involve neither Newton constant nor 
any other small parameter. In its turn, the local value of the scalar field C has a direct impact on the local values of 
the energy density and pressure as well as on the fermion mass. For this reason, first of all, we have to be sure that 
such ^ dependence of observables does not contradicts to experiments. We will make this in three steps. 

Step I. Let us start from the case of the absence of massive fermions. It follows from the constraint that in such a 
case C is a constant 

Then the effective cosmological constant in the fermion vacuum 



which is positive provided 

h > - 

sM4 



b>^.- (30) 



This inequality together with the one in Eg. (1^01) means that the integration constant sM^ > 0, that is s = 1. Since 
the measure $ is by definition sign indefinite, the scalar field ^ may be sign indefinite too. In some of the TMT 
models, this effect may be physically important [27| , [42j . But here we restrict ourself[43j proceeding in the region of 
the parameter space where Co > 0; that is 

Then (o ^ b and Eg. dTO]) may be completed by adding (q: 

Co-Ci-C2-b^k^h. (32) 

A tiny value of the vacuum energy density p^ac ^ M^/Ab may be achieved by the choice of a huge value of the 
dimensionless parameter b. For example, if M ^ (lO'^GeF)^ then for p^ac to be of the order of the present day vacuum 
energy density, the parameter b must be of the order b ~ 10^". One should stress that the choice of a huge value of 
the parameter b has nothing common with the fine tuning in its usual understanding in field theories (for details see 
item IV in the Discussion section). 
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Using Eas. ([28|) . pO)) and (ITTl) one can rewrite the constraint ([27)1 in the form 



^ (Co C) _ a(/) . a(/) „ p(noncan) _ [can] /„„n 

where the canonical fermion energy density reads 

p'r'-HOn- -(C)=A^(^^|(^^p5)T7i; /^>0- (34) 

and for the sake of simphcity we write down the contribution of a single fermion alone. 

It is very important that due to Eq. (j32p based on our basic assumption concerning the parameters of the model, 
Eq.®, the l.h.s. of the constraint has generically the same order of magnitude as Pvac, Eq. p5)) . if no a special fine 
tuning of the parameters is assumed. But clearly the l.h.s. of the constraint equals zero in the vacuum. 

Step II. Let us now check that the C, dependence of the fermion masses, Ea. ([M|) . has no observable effects in 
laboratory conditions. To show this, it is enough to notice that in all known so far particle ( as well as nuclear, 
atomic, etc.) physics measurements we deal with interactions which yield such localization of fermions that the local 
fermion energy density is many tens orders of magnitude larger than the vacuum energy density. In what follows, for 
short, we will refer to such a case simply as the high fermion density case. Hence if we want to estimate how the 
C dependence of the fermion masses can affect the values of observable mass one should compare m(C) in two cases 
with different but high fermion densities. 

In order to provide the validity of the constraint in the case of high fermion density of a nonrelativistic fermion 
one should find an appropriate value of C as the solution of the constraint. With the conclusion that the l.h.s. of the 
constraint is generically < p^jac and with knowledge of the obtained above expression (j29p for p^aa it is evident that 
the structure of the functions involved in the constraint does not allow to solve the problem by choosing |C| to be 
much bigger or much smaller than h. The only possible way to realize the balance between two sides of the constraint 
is to allow for C, to be very close either to Ci or to ^2- Then the factor Z{C,) in the r.h.s. becomes very small and it 
is able to compensate the large value of p^'^'^"-' such that the l.h.s. of the constraint remains ~ Pvac, as we discussed 
at the end of the step I. So, let us choose C, to be very close to Ci and represent it in the form Q = C,i{\ — 5) , where 
5 <C 1. Then the constraint may be written as following: 

M^(C0 - Cl) ^ '^Cl(Cl-C2) (can) . . 

(Ci+6)2 2(Ci + fc)(Ci+/i)^^ ^ ' 

For our purpose it is enough to know only the order of magnitude of 5. On the level of the orders of magnitudes, 
it follows immediately from Eq. ([55| that 

S ^ Pvac/ P'^/"''^ ■ (36) 

If we compare two cases of different high fermion densities P^(j^^^ and p\j^2) '^i^h the appropriate values of Si and 62 
then the relative change of the fermion masses is of the order 

\^^16,^62\<^, (37) 

Pf 

which is obviously unobservable. 

Step III. After studying two limiting cases of fermion vacuum and high fermion density, we are able to the search 
of the physically permitted interval of C,. In the case of cold fermions, it follows from Eqs. ((T^ and that the total 
energy-momentum tensor may be written in the perfect fluid form 

T^" = {P + P)u^'v!' - pg^'' = p^'^'^u'^u'^ + ktot'g^" (38) 

where is the four- velocity of the cold fermion and the following expressions for the total energy density and pressure 
read: 

(can 



p = pp^-f Atot, (39) 



p = -Ktot{0- 



(40) 
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The remarkable feature of TMT is that by virtue of the constraint (|33l) , the total energy density as well as the total 
pressure in the presence of cold fermions may be expressed in the form where they are functions of ^ alone: 



Af4 

P = 



P - 



lib-Co)+(: 



(41) 



(42) 



Our analysis is based on the following natural assumption: There should be possible transitions from the high 
density of cold fermions n — vj/'vl/' (when ( is very close to the value ^ or C2) to the low fermion density ended with 
asymptotic transition p — > p^ac as C — >■ Co- We suppose that C, as the solution of the constraint ([33]), is a continuous 
function of n evolving from C ~ Ci (or C ~ C2) to the regime C — > Co- In the course of the monotonic decay of n, the 
total energy density p must be continuous and positive function of C- This means that C cannot cross over values 
Ci and C2J where Z{() equals zero and p is singular. As it was already supposed by Ea. ((TT|) . C cannot cross also 
over the values —b and —k because the transformation to the Einstein frame, Eq.([S]), as well as the fermion mass, 
Ea. p4p . become singular and the constraint turns out to be senseless (it looks then as an equality of finite and infinite 
quantities). 

An additional information one can obtain from the constraint ()33|) taking into account that the fermion number 
density 



\ (C + &)^(C-Ci)(C-C2)- ^ 

is positive. Hence the physically permitted intervals for C must satisfy the condition: 

(C-Ci)(C-C2)(Co-C) >0. (44) 

Notice that by definition C2 < Ci- Now it is easy to check all possible intervals for C, and it turns out that there are 
only two regimes for evolution of C which satisfy to our natural assumption 

The case (A): C2 < Ci < C < Co; (45) 



The case (B): C2 < Co < C < Ci; (46) 

and the broad region in the parameter space exists where such regimes are possible together with the condition p > 0. 
In the case (A), the noncanonical fermion pressure 

pinoncan) ^ ^j^.j^ ^^^^ pinoncan) ^ ^ ^^^^ 

cases, in the course of the afore-mentioned monotonic decay of n, C remains of the order of the parameters b ^ k ^ h 
which means that the l.h.s. of the constraint (1331) has generically the same order (or very close to) as pvac, and hence 

,1 . T 1 r r^(noncan) i (noncan) 

the same is valid for Fj: ' and : 



I i (/) I I Yj(noncan) \ \ [noncan 



I [noncan), , i'A'r\ 
\Pf \ ^ -l^tot ^ Pvac (47) 



Only in the fermion vacuum A^JJ^ — and Atot — Pvac- Such the narrow intervals for values of C, A^Jyii ^tot 
explain why the noncanonical fermion energy density and pressure are unobservable under regular physics conditions, 
but in the cosmology of the late time universe they become important. This result follows from the constraint where 
the function Z{() plays the role of "self-locking retainer" permitting the narrow intervals for values of the listed 
quantities. 

V. COSMOLOGICAL AVERAGING IN THE LATE TIME UNIVERSE AND THE COSMOLOGICAL 

MODEL 



Let us now return to the end of Sec. Ill where we discussed the way to describe the evolution of the local energy 
density and pressure as governed by the geometrical, nondynamical scalar field Ci^)- In turn, the latter is evaluated 
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by numerical solution of the constraint (j27p as a function of the local magnitude of the fermion density. When trying 
to formulate a cosmological TMT model we have to answer the question whether the constraint, being very nonlinear 
in C algebraic equation, preserves its structure after cosmological averaging, that is it also has the form as in Eq. (|33|) 



(Co C) j=,{noncan) f){noncan) . (/) , 

{( + by ~^ ' ^ =^dyn^ (48) 

where now averaged value of pj"°"^°") ■ ^ jg cosmological averaged of its local space-time values 

C,{x) and now ^ is a function only of the cosmic time. 

Solution to this averaging problem turns out to be simple if we take into account the observed inhomogeneity of the 
universe: the existence of regions with clumped matter and domains where matter is very diluted. For example, in 
the large scale structure this is manifested in the existence of filaments and voids. The characteristic feature of such 
a structure is that, on each level of scales, the volume of the low density domains is tens orders of magnitude larger 
than the volume of regions with clumped matter. On the other hand, due to the self-locking retainer effect (see the 
end of Sec. IV) the relative differences both in values of C, and in values of pj"""^"") jj^ these two types of regions do 
not exceed one order of magnitude. Hence the cosmological average values C, and pj""""^"") -^rit]^ very high accuracy 
coincide with the values of <^ and pj""""^"") jjj the maximal volume domains of low fermion density. Therefore the 

averaged constraint indeed has the form of Eq. (j48p . where C and p("°"'^''") practically equal to the appropriate values 
in the maximal volume domains of low fermion density. 

In the cold matter dominated epoch, in the course of the cosmological expansion the universe enters the stage 
when in the maximal volume domains of low fermion density, p^^"'"'^ becomes less then Atot (C) • Then the regions of 
the universe populated with the maximal volume domains of low fermion density start to expand with acceleration 
which is, properly speaking, the transition from the cold matter dominated epoch to the DE dominated epoch. We 
come to the quite natural conclusion that the accelerated cosmological expansion of the universe can proceed due 
to accelerated expansion of voids and other low matter density domains. It is evident that such a manner of the 
accelerated expansion becomes more and more dominant in the course of evolution of the late time universe. Some 
speculations concerning possible approach to resolution of the puzzle why the galaxy formation epoch was about the 
transition to the DE domination epoch are presented in the item VII of the Discussion section. Here we want only 
to note that if the described mechanism of the accelerated expansion is right, it means that the transition to it could 
happen only after the voids are formed. But the latter must be accompanied by simultaneous formation of clustered 
matter. 

In the present paper we do not address the averaging problem of Einstein equations 35]. Since the Einstein Eas. (jl2p 
has the standard GR form, the first Friedmann equation in the spatially flat FLRW universe is as usual 

zp^b/' '''' 

where p should be the averaged total energy density. 

We would like to apply our model to the scenario which starts from the cold matter (dust) domination era, that is 
when the (dark) matter energy density is larger than the DE density. Let us suppose that the dark matter consists of 
fermions. Regarding the galaxy dark mater as the cold one we assume that its local conventional (thermal) pressure 
is negligible with respect to its local energy density. Then the local noncanonical energy density and pressure of 
the clumped fermionic dark matter being of the order of p^ac (see Eq. (|47|) ') are certainly negligible as well. It was 
shown in the step II of Section HI that the C dependence of fermion masses becomes negligible in space regions where 
fermion energy density is much larger than the vacuum energy density, and values of C are there very close to Ci or C2 • 
Galaxies and even filaments can be approximately regarded as such regions. Then the cosmological averaged energy 
density of the clumped matter scales as and we will describe it in the form 

Pel = -f (50) 

where M^j is the phenomenological parameter [i^l which allows to fit the duration of the cold matter domination era. 
Averaging of the second equation in ([55)) yields 



A^f„=^(c)p-r^ (51) 
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where the cosmological averaged of the canonical energy density of the cold fermion matter disposed in the maximal 
volume domains of low fermion density are obtained by averaging Ea. p4p : 

where n = n^/a^. Therefore p^"^""-* scales in a way different from \/a^. 

The local behavior of the the total variable CC term in the presence of fermions Tlt\x), Egs-dH]) and which 
as we hope is able to mimic the effect of DE, is governed by In the FLRW universe, the averaged components 

of T^it^ are obviously governed by C, and they read 

P(K) = -P(A) = -Atot(C), (53) 

where 



a,<,,(c) = a,,,(c)-a(^)„= ^'^ 



(54) 



and Ae// is determined by Eq. p3| . 

We would like to stress here that solving the constraint (|48| for C we obtain it as a function of n: C, = C(^)- We 
come to the crucial role result of the present paper: Atot is a function of the number density n of the cold fermion 
matter disposed in the maximal volume domains of low fermion density: 

Atot = Atot (n) (55) 

and the dynamical DE effect describing by Atot exists with no DE degrees of freedom. 
Ignoring the clumped matter pressure, we have that the total averaged pressure 

P = P(A) = -Atot- (56) 

The total (canonical) matter energy density consists of contributions of the clumped (dark) matter and the canonical 
energy density P(f^can) of the cold fermions disposed in the maximal volume domains: 

-{can) lrr,\ 
Pm = Pel + Pf (57) 

The total averaged energy density, which in the flat universe is the critical one, is then 

P = Pin+P(K), (58) 

and fractions of the normal matter energy density and of the one which mimics the DE density are then defined as 
usual: 

an = — and nDE = ^ (59) 
P P 

Making use of the constraint. Eg. (1481) . one can rewrite p and p in the form where the contributions from the maximal 
volume domains depend only on the averaged scalar field C,: 



Pel 



2^ - Z(C) 



(60) 



(61) 



The constraint (j48|) is the 5-th degree algebraic equation which determines C as an implicit function of pn. Since 
the constraint is valid for any value of t, the time derivatives of both sides of the constraint are also equal at all t. By 
virtue of this evident fact and using the Friedmann equation (j49p we obtain the first order differential equation for C 

(62) 



Bt-W,iCo-C)V^P 
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where 

+ + + (63, 

It is much more convenient to work with the differential equation than to proceed with the constraint itself. 
Notice that 

p + p = Pm, (64) 

where p,„ is determined by Eq. (j57|) . Then the total averaged energy- momentum tensor may be written in the perfect 
fluid form 

T'^'^ = p^u'^u"^ + Atotr" (65) 

where = diag{l,—a?,—a?',—a'^), the four-velocity of the cold fermion gas in the co-moving frame is — 
(1,0,0,0) and Ktot is determined by Eq. ([54)) . The structure of the cosmological model under consideration is very 
unusual. For example, the fermion masses and Ktot are C, dependent and evolve in the course of cosmological expansion. 
In spite of this, the standard GR equation of the energy-momentum conservation 

p=~?,-{p + p) (66) 
a 

is satisfied. This statement may be checked by means of bulky enough calculations making use p and p, Eqs. (|Bn| . 
(|6T|) . the constraint (|48l) and Eas. ([62l) and (|49l) . Suppose for the moment that we have solved the system of Eas. (|49)) 
and (|62|) . Then we have p and p as functions of C which in turn is a function of cosmic time (or of the scale factor). 
These functions may be treated as the function p = p{p) given in the parametric form. We would like here to stress 
again that even though we deal effectively with a kind of barotropic fluid, the model does not involve any degree of 
freedom intended to be dynamical DE. 

To solve the system of Egs. p^ and (p^ we need an additional information which allows to choose the initial 
condition for C,. This information one can obtain from the constraint, Eq. (|48|) . taking into account that at the very 
beginning of the scenario described by our model, i.e. in the cold fermionic matter domination era, m{Q)n is many 
orders of magnitude larger than the l.h.s. of the constraint p8)) which is of the order of (or very close to) the present 
vacuum energy density (see Eq. (|47l) '). This is possible if Z[C,), Eq.(fT7|). is very close to zero in the cold matter 
domination era. Taking into account the physically permitted intervals of C, Eqs. (H5t and (H^ . we must choose the 
initial value Ci„ to be very close to Ci. Then the EoS w = p/p is very close to zero and negative (see Eqs. (|5(I| . (pT|) 
and (HZ])). In the opposite limiting case, when a{t) — >■ oo, it follows from the constraint that — > Co, where Co is 
the vacuum value of C, defined by Eq. ipS)) . One can immediately see from Eqs.(|Bn| and (|5T|) that w asymptotically 
approaches w = —1. 

Using the results of the qualitative analysis in Sec. IV (Step III) one can classify the dynamics of the late time 
universe in the correspondence with the region of the parameters space: 

• The LA type of scenario with w > —1, when —b < Ci, —k < Ci and —h < Ci < Co and C satisfies the inequalities 

(A) , Eq.m, 

C2 < Ci < C < Co; (67) 

• The LB type of scenario with w > ~1, when — 6 < Co, — fc < Co and — /i < Co < Ci and C satisfies the inequalities 

(B) , Eq.®, 

C2<Co<C<Ci; (68) 

The LA and LB types of scenarios are very similar: the total energy density p monotonically decreases to Pvaa 
Ea. (P^ . and the EoS monotonically decreases to w — —1. The difference consists in the opposite behavior of 
C: in the LA type of scenario C increases approaching the vacuum value Coi while in the LB type of scenario 
C decreases approaching Co- Such behavior is dictated by Eq. (p^ : in the LA case C > 0, while in the LB case 

C<o. 
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• The II, Phantom-like type of scenario, when —b < Coi ~k < Co: Co < ~h < Ci and C - as in inequalities (j68l) 
It follows from Eg . ((62)) that ^ < 0, i.e. C, decreases approaching Co- It is easy to see from Eas. ([52)) . p8|) . ((5T 
and dni) that 



sign (to(C)'t.) = sign 



-(Co-C) 



sign{C + h) (69) 



In the course of decreasing, C crosses the value —h. At this moment, both the effective fermion mass and 

Eg. ([5^ . cross zero values and become negative, but pj^""-* ~^ P{A) remains positive provided that the parameters 
are chosen such that, in addition to the listed conditions, 

i(?^-Co) + C+^>0 (70) 

not only for —h < C < Ci but also when Co < C < By means of numerical estimations one can check that 
there is a broad range of parameters where all the above conditions are satisfied. 
Taking into account that + P(A) = we have for the EoS 

w+i^kpci+p'^rh (71) 
p 

At the moment when ^ = —h, the fermions possess only the noncanonical pressure and energy density and 
their EoS coincides with the cosmological constant EoS. Such a vacuum- like fermion state has nothing common 
with usual conception of what is a particle. One can suppose that from all known particles, only neutrinos in 
maximal volume domains of low fermion density are able to reach such exotic state. This state may be treated 
as the limiting realization of the idea of the neutrino dark energy[l3i-[i3|- If the cosmological averaged energy 
density of the clumped matter pd is less than the DE density p^^) (that is in the agreement with the present 
day cosmological data) then soon after C crosses the value —h, w should cross the value w = —1. Further 
cosmological expansion proceeds with |m(C)fi| (that causes C — Co —J' 0) and pd — >■ 0. Therefore p should 
pass through minimum that is followed by monotonic increasing with p — ^ p^ac from below and w — >■ — 1 from 
below as well. Notice that here pyac is determined by Eg. d^ and, as we have seen in Sec. IV, Eg. (|T7|) . it is of 
the order of (or very close to) the present vacuum energy density. Such an asymptotic behavior belongs to the 
pseudo-rip type of the future of the universe, according to the classification of Refs.fSG],!^. Such a phantom-like 
cosmology of the very late time universe is realized here without phantom scalar field (15i] as well as without 
modification of the Friedmann equation[l^ based on the braneworld model[3|. However since the cosmological 
dynamics of the model is governed by the effective barotropic fluid (see the discussion after Eg. lpS)) ) all problems 
of singularity and instabilities (with respect to cosmological perturbations) related to the possibility of crossing 
the phantom divide [3^,[3^ are relevant here and have to be studied. This will be done in a separate publication. 

• The III, Phantom-like scenario with unlimited grows of |m(C)| is realized with the same conditions as in the 
II type of scenarios except for the very special fine tuned choice Co = In this case, the scenario is similar 
to the II type of scenarios but one can see that in the last stage, as a — ?> oo and C ~^ Co = ~k, the neutrino 
mass is negative and behaves as m(C) ^ (C + fc)~^; the fermion number density n, Eg. (|43)) . behaves as 
n ~ 1/a^ ~ (C -I- k)^. Therefore the canonical neutrino energy density. Eg. ([5^ . being negative approaches zero 
as pj"""'-' ~ (C + k)^ ^ while p(A) ~^ Pvac- Notice that this qualitative analysis is consistent with Eg. ([5^ 

since in the regime when C being larger than —k appears to be also close to —k then the equation yields C — 0~, 
i.e. C is forced to tend to ~k. 

VI. THE MODEL RESULTS COMPUTED NUMERICALLY 

The results of numerical computations are presented in Figs. 1-20 where we have used the dimensionless units 
obtained by the following changes: 

-^C; -^fc; -^h- ^^Ao; -^Mti^Mtu 
b ^ b b bM^ ' M4 ''^ 

all others energy densities and pressures: -r^rP — > p and -rrjP — > P- (72) 
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• The typical features of the LA type of the cosmological evolution scenario are illustrated in Figs. 1-7 with the 
choice ofjthe parameters Ap = -1.5, k = 2, h = 0.3, = 10"^. Then Ci ^ 2.3114, (2 = -1.2114. The initial 
value of C is chosen to be (in = 2.312. 

• The typical features of the LB type of the cosmological evolution scenario are illustrated in Figs. 8-13 with the 
choice of the parameters: Aq = -0.2; k = -0.3; h = -1; = 5 • 10^^. Then Ci = 3.0595; C2 = -0.3595. The 
initial value of C is chosen to be (in = 3.0594. 

• The typical features of the II, phantom-like type of the cosmological evolution scenario are illustrated in Figs. 
14-19 with the choice of the parameters: Ao = -0.2; k = -0.3; h = -6; = 5 • 10"^; Then Ci = 18.2267; 
(2 = —0.5267. The initial value of C is chosen to be Ci„ = 18.2266. The effect of the change of the neutrino 
mass (in the maximal volume domains of low density) that yields the phantom-like cosmological dynamics is 
iUustrated in Fig. 17. 

• The typical features of the III, phantom-like type of the cosmological evolution scenario are similar to those of 
the II type. The only difference consists in the fine tuning of the parameters such that Co = which yields 
the unlimited approach to —00 of the neutrino mass in the maximal volume domains of low density. This effect 
is illustrated in Fig. 20. 

The results of the numerical solutions confirm the classification and properties of possible cosmological scenarios 

obtained by the qualitative analysis in end of Sec.V. 

The key point is the confirmation of the self- locking retainer effect: as the total energy density decays from a very 

large value at the cold matter dominated epoch up to the very low value at the DE dominated epoch, the changes 

- ...(f) 
of the averaged scalar C and the averaged value of the dynamical fermionic A term A^^^^ do not exceed one order of 

magnitude. Due to this effect, the interval of the change of the total variable CC Atot during the studied history of 

the universe does not exceed one order of magnitude too. Remind that the effect of the dynamical DE described by 

Atot is generated and driven by the cold fermions whose density changes from the value at the cold matter dominated 

epoch up to the density of cold neutrinos in voids which tends to zero with the cosmological expansion. 

Another interesting result of the numerical solutions is that in all possible scenarios the transition from the cold 

matter dominated epoch to the DE dominated epoch occurs rapidly about the red shift z ~ 1, and this is realized 

without any special choice of the model parameters, integration constant and initial conditions intended to achieve 

this result. 



VII. SUMMARIZING THE RESULTS AND DISCUSSION 

/. Fermions in TMT. One of the peculiarities of the TMT consists in the perfectly novel features of massive fermions, 
and the geometrical scalar field (, Eq.®, plays the main role in this respect: 

1. The fermion mass is C dependent, while the local value Ci^) is determined via the constraint as a function of the 
local fermion density. As the local fermion energy density p/ is many orders of magnitude larger than the vacuum 
energy density Pvac, the fermion mass is constant up to corrections of the order of Pvac/pf, and the local value C{x) 
must be very close either to Ci or to C2 (see the step II of Sec. IV). Thus, such "high density" fermion states belong to 
the regular particle physics situation. 

2. As p/ becomes comparable with p^ac, the fermion mass becomes C dependent, and due to the constraint this 
means that the fermion mass depends upon the fermion density. For short, we called such kind of fermion state as 
the low density one. 

It follows from the last two items that in order to describe the fermion state in the quantum field theory formulated 

in the framework of TMT, in addition to the usual quantum numbers one should add the value of C- 

...(f) 

3. In addition to the canonical energy-momentum tensor, fermions possess also the dynamical fermionic A^^^^ term. 

Eg. (1211) . and therefore the appropriate noncanonical pressure P(f^noncan) = -^iyL and the noncanonical energy density 
P(f,noncan) Satisfy the equation P(f^noncan) = -P{f,noncan)- Wc havc shown that p^f^noncan) and p^ac are always of the 
same or very close orders of magnitude. For that reason the effect of the noncanonical energy-momentum tensor is 
unobservable in all particle physics measurements till now but it has a crucial role in the cosmology of the late time 
universe. 

//. Dynamical DE imitated by the geometrical effect of TMT in the presence of massive fermions. The matter 
content of the 4D TMT model studied in the present paper is reduced to massive fermions though radiation could be 
added without altering the results because gauge fields do not contribute to the constraint. The important feature of 
the model is that it does not involve any fiuid intended to describe the DE. In the case of the absence of fermions, it 
follows from the constraint that C = const and then the model contains the constant CC. However, in the presence of 
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fermions, on account of the constraint, ( becomes a function of the local fermion density. In its turn, the constant CC 
turns into the C dependent CC, i.e. on account of the constraint, the CC becomes dependent of the fermion density. 
This CC together with A)j^'„ form the effective variable CC, Atot, Eq.([24| depending upon the fermion density. Note 

that since Atot contains the term AJy^^ it would be strange, as usually, to treat p(A) = Atot as the variable vacuum 
energy density. It seems to be more natural to regard Atot as the effective dynamical DE imitated by the geometrical 
effect of TMT in the presence of massive fermions and governed by the fermion density. 

///. Cosmological averaging and cosmology of the late universe. We have concentrated on the cosmological outputs 
of the key effect of the model: generating of the variable CC Atot which: (firstly) mimics the effect of DE and 
(secondly) evolves as a function of the cold fermion density n. We have shown that the constraint keeps its form after 
the cosmological averaging and in the course of the cosmological expansion, the cosmological average of Atot with the 
very high accuracy may be regarded as a function of the value of n in the maximal volume domains with the low 
fermion density. In other words, one can claim that the effective cosmological dynamics imitating the DE effect is 
governed by the cold fermions density in domains like voids and supervoids. We have proposed to associate these 
fermions with neutrinos because it seems to be natural to assume that, in the course of cosmological expansion, from 
all known fermions only nonrelativistic neutrinos are able to survive without interactions up to the state when their 
canonical (thermal) pressure is negligible in comparison with their noncanonical pressure. If this assumption is true 
then the DE effect is driven by neutrinos and the neutrino DE is the essential fraction of the total effect. 

We have find the permitted by TMT intervals of ^ where cold fermions can evolve continuously from the high 
density to the low fermion density state ended with asymptotic transition p Pvac- Studying cosmological evolution 
of the late time universe in these intervals we have shown that there exists a broad range of parameters where the 
rapid transition (at z '-^ 1) from the cold matter domination era with EoS w w to the EoS w « —1 occurs. In this 
paper we did not fit parameters for better concordance with observational data, but instead of this we concentrated 
on the field theory argumentation in favor of the geometrical effect imitating dynamical DE. In particular we did not 
address to the cosmic coincidence problem ("why now?"), and the clustered dark matter density has been accounted 
phenomenologically. 

IV. The value of the vacuum energy density The effect of DE was realized here on the basis of the first principles 
of TMT. The only additional quantitative assumption is that the dimensionless parameters 6, k and h have the same 
order of magnitude. As we mentioned in the step I of Sec. IV, the observed tiny value of p^ac is achieved if these 
parameters are huge numbers. We want to argue here that the latter choice is very much different from what is 
usually called a fine tuning in resolution of the old or new CC problem [ioj. When the CC problem one hands over 
to particle field theory, the generally accepted understanding of the required subtraction or cancellation mechanism 
consists in a fine tuning of coupling constants and masses in similar terms in the Lagrangian, like e.g. as it is done 
in SUSY and SUGRA. The way large numbers we would like to appear in the TMT action is absolutely different. 
Inserting the factor h in front of S in Eq.Q we want to stress that the huge values of 5, k and h can be equivalently 
treated as an extremely small value ^ 1/6 of the coupling constant of the (gravitational + matter) Lagrangian Li 
to the measure see also Eq.([T|). Thus the tiny parameter 1/5 has no relation to any hierarchy or fine tuning in 
the particle field theory but rather it may indicate that the emergence of the Lagrangian to the measure $ coupling 
in the action results from breaking of some unknown symmetry as yet. One can speculate, for instance, that it is a 
quantum gravity effect [45]. _ 

V. Generalization to the Standard Model. In our earlier papers we have shown that the gauge fields may be 
added to the model without changing the constraint. The nonabelian symmetries are also compatible with TMT. 
Therefore the matter content of the model may be extended up to that of the standard model. Note that all the TMT 
effects emerging in the presence of fermions become unobservable in the regular particle physics conditions. We would 
like to stress again that results of the present paper allow us to claim that the effect of DE may be realized with the 
only matter content of the standard model and without any additional dynamical degrees of freedom intended to get 
this effect. 

VI. Phantom-like cosmological scenario. In the broad region of the parameter space, the model admits a phantom- 
like cosmological scenario where the crossing of the phantom divide is caused by changing the sign of the effective 
neutrino mass in the course of the cosmological expansion. As the result, at this moment the neutrinos energy density 
and pressure satisfy exactly the CC EoS that means disappearance of neutrinos as particles (excitations) in a sense 
this term has in the quantum field theory. Such a limiting realization of the idea of the neutrino DE deserves more 
detailed investigation. Even though the sign of fermion mass in the relativistic quantum mechanics is a question of 
choice, it remains unclear the physical meaning of the new cosmological dynamics effect we observe in our model: 
"degeneration" of neutrino with positive effective mass with subsequent "regeneration" with negative effective mass. 
Nevertheless this strange effect, from our viewpoint, seems to be less dangerous than the use of the phantom scalar 
field [i^. We should note also that in this paper we restrict ourself with studying massive fermions without taking 
into account the Higgs mechanism. Preliminary study shows that in TMT, in the context of cosmology, a possible 
a back reaction of fermions on the vacuum expectation value of the Higgs field may be relevant. In this respect the 
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model of the present paper may be regarded as a toy one. It would be interesting to investigate in the future what 
kind of effect such a back reaction may have on the change of the sign of the neutrino mass and on the possibility of 
the phantom-like scenario. 

VII. Fermion dark matter candidate? One can finally propose the hypothesis concerning the nature of the dark 
matter. Creation of fcrmions in regular particle physics conditions implies that the fermion emerges in the "high 
density" regime. As we have mentioned in the item I of this section, in such a case the fermion may be in the state 
with very close either to Ci or to (^2- In our earlier papers [26] concerning studying the Standard Model in the 
framework of TMT, we explored the idea that fermion states with C very close to Ci and C2 might be associated with 
the first and second fermion generations respectively. However in the present paper we run into a drastic difference 
in features of the ^1 and ^2 fermion states. In fact, we have elucidated that if fermions have been created in the state 
with ^ very close to then they are able to evolve to the arbitrary low energy density. If however fermions have 
been created in the state with ^ very close to C2 then it follows from Eq.(|33|) that only in two following cases it is 
possible: (1) Co < C2 < C < Ci and (2) — & < C < C2 < Co < Ci (the case (2) involves an additional requirement on 
the parameter b which is satisfied in the numerical solutions of Sec. VI for the scenarios lb, II and III). One can check 
using Eq.(|62|) that C is forced to oscillate in the interval C2 < C < Ci in the case (1) and in the interval — 6 < C < C2 in 
the case (2). Such regimes should perhaps be ended at some intermediate respective values, but for us it is important 
here that the fermion starting from such a state is unable to evolve to a state with arbitrary low energy density. 
We come to this conclusion in the FLRW metric but one can expect a similar result e.g in a spherically symmetric 
gravity problem. One can also expect that there should be a mechanism which prevents a possibility for the fermion 
in the states of the cases (1) and (2) to be in the same space-time region as well as for each of them to be in the 
same space-time region with another fermion whose state starts from C very close to Ci- Suppose that the dark 
matter consists of the same fermions (from the viewpoint of the quantum field theory) as the visible fermion matter. 
Our hypothesis may be formulated in the following form: the fermion in the state with C very close to Ci is the 
visible fermion while the same fermion but in the state as in the cases (1) or (2) is a candidate for the dark matter 
particle [13] • If this hypothesis is true then there should be a new mechanism responsible for spatial separation of 
visible and dark matter, that could highlight the structure formation in a new fashion. It is worth to recall that 
analyzing the procedure of the cosmological averaging we came to the conclusion that with the appearance of voids 
and other low matter density domains, the averaged value C with very high accuracy becomes equal to the value of C 
in such domains. But low matter density domains constitute the part of the universe where Atot may be larger than 
the canonical fermion (neutrino) energy density. In such a case accelerated expansion of these domains produce the 
effect of the accelerated expansion of the visible part of the universe. Since emergence of such domains is apparently 
accompanied with formation of regions with clustered matter one can hope that this is a way to explain the puzzle of 
the approximate time coincidence of the galaxy formation epoch with the transition to the DE domination epoch. 
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FIG. 1: In figures 1-7 the LA type of scenario of the cosmological evolution with parameters: Ao = —1.5, k — 2, h — 0.3, 
M^i = 10"^. Then Ci = 2.3114, (2 = -1.2114. The initial value of C is chosen to be Cm = 2.312. In this figure: total EoS. 




-8 -7 -6 -5 -4 -3 -2 -1 

log a 



FIG. 2: The LA type of scenario. Total energy density. 
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FIG. 3: The I. A type of scenario. The averaged dynamical fermionic A term, Eg 1541 vs (. The result of this numerical solution 
confirms the analytic estimations made in the end of Sec. IV and formulated under the name the "self-locking retainer" effect. 




FIG. 4: The LA type of scenario. The averaged dynamical fermionic A term, Eg 1541 vs In a. 
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FIG. 5: The LA type of scenario. Mass of the cold neutrinos in the maximal volume domains of low density: m vs Ina and m 
vs C 
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FIG. 6: The LA type of scenario. Evolution of the averaged value f that with very high accuracy coincides with the value of 
in the maximal volume domains of low fermion density. The result of this numerical solution confirms the analytic estimations 
made in the end of Sec. IV and formulated under the name the "self-locking retainer" effect: as the total energy density decays 
from the value ~ 10^'' at the cold matter dominated epoch with ~ up to the value ~ 10~^ at the DE dominated epoch 
with w ~ —1 (see Fig. 2), ( changes only from ( ~ 2.3 to the fermion vacuum value ^ — (o ~ 4: defined by Ea. (|28p . 
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FIG. 7: The LA type of scenario. f2 vs Ina where fractions of matter Qm (the black dash line) and effective DE ^Ide (the blue 
solid line). 
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FIG. 8: In figures 8-I4 the LB type of scenario of the cosmological evolution with parameters: Ao = —0.2; k = —0.3; h = 
M^i = 5 ■ 10"^ Then = 3.0595; C2 = -0.3595. The initial value of ( is chosen to be C„i = 3.0594. In this figure: total EoS. 
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FIG. 9: The LB type of scenario. The averaged dynamical fermionic A term, Eg 1541 vs ( and vs In a. The result of this numerical 
solution confirms the analytic estimations made in the end of Sec. IV and formulated under the name the "self- locking retainer" 
effect. 
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FIG. 10: The LB type of scenario. Total energy density. 




FIG. 11: The LB type of scenario. Mass of the cold neutrinos in the majdmal volume domains of low density: m vs Ina and 
m vs C,. 
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FIG. 12: The LB type of scenario. Evolution of the averaged value ^ that with very high accuracy coincides with the value of ( 
in the maximal volume domains of low fermion density. The result of this numerical solution confirms the analytic estimations 
made in the end of Sec. IV and formulated under the name the " self- locking retainer" effect: as the total energy density decays 
from the value ~ 10^® at the cold matter dominated epoch with w « up to the value ~ 10~^ at the DE dominated epoch 
with w ~ —1 (see Fig. 10), ( changes only from f ~ 3.06 to the fermion vacuum value ( — (^o = 1.4 defined by Eq. (|28|) . 
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FIG. 13: The LB type of scenario. 57 vs Ina where fractions of matter f2,„ (the black dash line) and effective DE Qns (the 
blue solid line). 
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FIG. 14: In figures 14-19 the II, phantom-Uke type of scenano of the cosmological evolution with parameters; Ao = —0.2; 
k = -0.3; h = -6; M^i = 5 • 10"*; Then Ci = 18.2267; C2 = -0.5267. The initial value of C is chosen to be On = 18.2266. In 
this figure: the total EoS typical for the phantom-like scenario of the late universe 
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FIG. 15: The II, phantom-like type of scenario. The averaged dynamical fermionic A term, Ea l54l vs C, and vs In a. The result 
of this numerical solution confirms the analytic estimations made in the end of Sec. IV and formulated under the name the 
"self-locking retainer" effect. 
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FIG. 16: The II, phantom-like type of scenario. The result of the numerical solution for the total average energy density 
p confirms the analytic estimations made in Sec.V: p passes through minimum that is followed by monotonic increasing with 
p — ^ pvac from below. Such a behavior of p is typical for the phantom like scenario where the very late universe evolves as in 
the pseudo-rip mode/ [36|. 
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FIG. 17: The II, phantom-like type of scenario. Mass of the cold neutrinos in the maximal volume domains of low density: m 
vs In a and m vs i^. It is seen that m changes sign and asymptotically approaches a finite negative value. 
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FIG. 18: The II, phantom-like type of scenario. Evolution of the averaged value f that with very high accuracy coincides 
with the value of in the maximal volume domains of low fermion density. The result of this numerical solution confirms the 
analytic estimations made in the end of Sec. IV and formulated under the name the " self- locking retainer" effect: as the total 
energy density decays from the value ~ 10^^ at the cold matter dominated epoch up to the value ~ 10~^ at the DE dominated 
epoch (see Fig. 16), ^ changes only from ~ 18 to the fermion vacuum value C = Co = 1-4 defined by Ea.(|28[l. 
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FIG. 19: The II, phantom-like type of scenario. vs Ina where fractions of matter {Qm, the black dash line) and effective DE 
(Qde, the blue solid line). The simultaneous effects of exceedance of the value 1 for Qde and lowering below the value —1 for 
flm result from changing sign of the neutrinos mass and their canonical energy density in the maximal volume domains of low 
fermion density. 
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FIG. 20: The HI, phantom-like type of scenario of the cosmological evolution with the fine tuned parameters such that 
Co = -k: Ao = -0.2; k = -1.4; h = -6; M*, = 5 • 10"*; Then Ci = 16.6481; C2 = -0.0481. The initial value of ( is chosen to 
be ("in ~ 16.647. In this figure we show the evolution of the mass of the cold neutrinos in the maximal volume domains of low 
density: m vs Ina and m vs ^. m changes sign from positive to negative and m — >■ —00. Cosmological evolution of all other 
quantities is similar to that of the II, phantom-like type of scenario presented in Figs. 14-16, 18 and 19. 
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